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Abstract 

A detailed study of critical spreading in the one-dimensional pair contact 
process is performed using a recently devised reweighting method. The re- 
sults confirm the validity of a generalized hyperscaling relation among the 
(nonuniversal) spreading exponents, and support the assertion that the criti- 
cal point location does not depend on the particle density <j). It appears that 
the exponents z and 5 + rj, once thought to be invariant, exhibit a systematic 
dependence on </>. 
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I. INTRODUCTION 



The model studied in this work is a not-to-distant cousin of directed percolation. In 
view of his extensive and continuing contributions to percolation theory, it gives me great 
pleasure to dedicate this work to my friend George Stell, who has taught me much about 
physics (and perhaps more about music). 

When viewed as a spatio-temporal process (the oriented axis representing time), directed 
percolation exhibits an absorbing-state phase transition, that is, a transition between an 
active state and one in which the dynamics is frozen. Critical phenomena at absorbing-state 
phase transitions are of longstanding interest in statistical physics [1,2], being found, for 
example, in models of epidemics [3], catalytic kinetics [4,5], surface growth [6], and self- 
organized criticality [7-10]. The best studied models of this kind, the contact process (CP) 
and directed percolation (DP), are relevant to experiments on interface pinning [11] and sand 
flow [12]. The transition between active and absorbing states arises out of a conflict between 
two opposing processes (e.g., creation and annihilation); when continuous (as is often the 
case), it falls generically in the DP universality class [13-15]. When two or more absorbing 
states exist and are connected by a symmetry operation, as in branching and annihilating 
random walks, the asymptotic dynamics is dominated by a domain-growth process, and a 
new kind of critical behavior appears [16-19]. 

Unusual critical behavior also appears in models that can become trapped in one of an in- 
finite number of absorbing configurations (IN AC). (More precisely, the number of absorbing 
configurations grows exponentially with the system size. There is no special symmetry link- 
ing the different absorbing configurations.) Models of this sort were introduced to describe 
surface catalysis [20-22]; their critical properties have been studied in detail by various work- 
ers [23-30]. In one dimension, the pair contact process (PCP) [24], and other models with 
INAC exhibit static critical behavior in the DP class [25,31], but the critical exponents 5, rj, 
and z, associated with the spread of activity from a localized seed are nonuniversal, varying 
continuously with the particle density (f> in the environment [25,26,30,32]. [These exponents 



are defined via the asymptotic (t — > oo) power laws: survival probability P(t) ~ t~ s , mean 
activity n(t) ~ t v , and mean-square distance from the seed R 2 (t) ~ t z ] This anomalous 
aspect of critical spreading for INAC can be traced to a long memory in the dynamics of 
the order parameter, p, arising from a coupling to an auxiliary field (the local particle den- 
sity, 0), that remains frozen in regions where p = [27,29,30]. Theoretical understanding of 
models with INAC remains incomplete. Muhoz et al. were able to construct a plausible field 
theory (a set of stochastic partial differential equations), for models with INAC, and to show 
that the static critical behavior is that of DP [27]. Formally eliminating the auxiliary field, 
they obtained a closed equation for the order parameter, in which a memory term appears; 
simulations of this theory reproduce the nonuniversal exponents observed in simulations of 
particle models [32]. On the other hand, the phenomenon of nonuniversal spreading expo- 
nents has so far resisted analysis via renormalization group or other theoretical methods. 

Studies of spreading in models with INAC have yielded several clues toward understand- 
ing nonuniversality [25,26]. First, the location of the critical point (the critical parameter 
value, p c in the PCP), does not depend on 0, even as the exponents vary. Second, the 
exponent z, and the sum 5 + rj, are roughly independent of 0. Third, the exponents 5, i], 
and z obey a generalized form of the hyperscaling relation derived by Grassberger and de la 
Torre for models with a unique absorbing configuration [26,33]: 

v + 5 + ^ = -. (1) 

u \\ 2 

For models with a unique absorbing configuration, the scaling relation 5 — f3/v\\ holds, and 
one recovers the original hyperscaling relation [34]. When the static critical behavior falls in 
the DP universality class, f3/v\\ = 6 dp independent of 0; in what follows we shall write Eq. 
(1) in the form rj + 5 + 5 DP = dz/2. Recently, Odor et al. reported a weak dependence of p c 
on the particle density in the environment [35]. In this work I describe extensive simulations 
of the PCP, using a recently-devised reweighting method [36], and improved data analysis. 
The results shed light on all three of the above-mentioned clues. 

The balance of this naner is devoted to defining' the model and simulation algorithm 



(Sec. II), some observations on the implications of a critical point shift for scaling (Sec. Ill), 
outlining the simulation method (Sec. IV), and presenting the simulation results (Sec. V). 
We close (in Sec. VI) with a discussion of our findings and their relation to previous work 
on nonuniversal spreading in models with an infinite number of absorbing configurations. 

II. MODEL AND SCALING PROPERTIES 

Jensen's pair contact process (PCP) [24], is an interacting particle system: a Markov 
process whose state space is a set of particle configurations on a lattice [37,38]. Here we 
consider the one-dimensional version: each site of Z d is either vacant or occupied. Each 
nearest-neighbor (NN) pair of occupied sites ( "particles" ) has a rate p of mutual annihilation, 
and a rate 1 — p of attempted creation. In a creation event involving particles at sites i and 
a new particle may appear (with equal likelihood) either at site i — 1 or at i+2, provided 
the chosen site is vacant. (Attempts to place a new particle at an occupied site fail.) In an 
annihilation event the sites occupied by a NN pair are simply vacated. The PCP exhibits 
an active phase for p < p c ; above this value the system falls into an absorbing configuration 
devoid of NN pairs, but that typically contains a substantial density, 0, of isolated particles. 
The most precise available estimate for the critical parameter isp c = 0.077090(5) [31]. (Here 
and in what follows, numbers in parentheses denote an uncertainty estimate in the last figure 
or figures.) Note that the active regime corresponds to p < p c . 

Previous studies leave little doubt that the static critical behavior of the PCP belongs 
to the universality class of directed percolation. Jensen and Dickman found that the critical 
exponents f3, 7, 6 (which govern, respectively, the stationary mean of, variance of, and initial 
decay of the order parameter), and v\\ and u± (which govern the divergence of the correlation 
time and correlation length as one approaches the critical point), are all consistent with DP 
values [25]. More recently, the order parameter moment ratios and cumulants were found 
to be the same as those of other models belonging to the DP universality class [31]. The 
process of critical spreading, i.e., the propagation of initially localized activity in a system at 



the critical point, is more subtle. Critical spreading in the PCP is studied by simulating the 
process at the critical point, with an initial condition generated by placing a single NN pair 
of particles (the "seed") into a configuration otherwise devoid of pairs. In this work, as in 
previous studies, the sites outside the seed bear a uniform particle density 0. (The particles 
may be placed in a random or a regular fashion.) The critical point for spreading, p c (<f>), is 
defined operationally as the value for which the survival probability, mean number of pairs, 
and mean-square distance from the seed follow power laws. All earlier studies, except that 
of Odor et al., have found p c (<f>) = p c , independent of 0. 

The surprising result of the critical spreading studies is that only for a particular value 
of 0, the so-called natural density 4> na t, do the spreading exponents 5, i] and z assume 
DP values [25]. The natural density is defined as the mean particle density in absorbing 
configurations generated by the process itself, at p c , starting from a homogeneous (e.g., fully 
occupied) configuration. One may, equivalently, define (j) nat as the particle density in the 
critical stationary state, in the thermodynamic limit. An environment with > <p nat should 
favor spreading (and one with < <p nat should hinder it), since the higher the particle 
density in the environment, the more pairs will be formed per creation event. 

A kind of spreading phenomenon also arises in the stationary state due to spontaneous 
fluctuations. In the critical stationary state, we can expect to find inactive regions of all 
sizes; the particle density in large inactive regions is <fi na t- When activity spreads into such 
regions, it should follow the same scaling behavior as critical spreading with nai . Since the 
exponents governing survival and growth of activity in the stationary state are subject to 
the scaling relations 5 = f3/v\\ and z = 2v±/v>\\, with f3, u± and u\\ taking DP values in the 
stationary state, it follows that the spreading exponents take their usual DP values as well, 
for = (j) na t- In an environment with ^ (fi na t, the advance of the active region is no longer 
equivalent to the spread of activity in the stationary state, and the spreading exponents 
are not constrained to take DP values. Of course, the interior of the active region must 
eventually relax to the critical stationary state, with a particle density (p na t, regardless of 
the exterior narticle densitv (h. 



In simulations one finds that for < (f) nat the exponent 5 is larger than the DP value, so 
that the survival probability decays more rapidly. r\ and z are smaller than the corresponding 
DP values, showing that diminishing the particle density hampers spreading, as expected. 
For > (fi nat these trends are reversed. The variations in 5 and r\ are quite dramatic: the 
former changes by more than a factor of four, and the latter by about two, as is varied 
between the extreme values of and 1/2. The changes in z and 5 + r\ are much smaller, 
amounting to about 5% over the full range of 0; it was suggested that the observed variations, 
which were not much larger than numerical uncertainties, could be attributed to corrections 
to scaling, suggesting the appealing simplification that z and 5 + i] are in fact independent of 
[26] . Since z and 5 + 77 have only to do with surviving trials (the latter governs the growth 
in the number of pairs in such trials: n s ~ t 5+v ), lack of dependence of these quantities on 
would imply that the latter only affects the survival probability, and not the asymptotic 
scaling of surviving trials. 

Two observations should be made before discussing further details. First, the above pat- 
tern of scaling properties (DP static behavior, nonuniversal exponents taking DP values only 
at the natural density, no apparent shift in p c , and near-constancy of z and 5 + r)), has been 
confirmed in other one-dimensional models with INAC (the so-called dimer reaction [25], 
the threshold transfer process [26], and an evolution model with real- valued site variables 
[39]), in models with a long memory of the initial configuration [29,30], and in a field the- 
ory intended to describe models with INAC [27,32]. Second, spreading in two-dimensional 
models with INAC appears to be more complex [29], and may not exhibit power- law scaling 
for all initial conditions [30], although the static critical behavior again falls in the DP class 
[28,40]. 



III. CRITICAL POINT SHIFTS AND SCALING 



How will shift in the critical point affect the spreading process? Before attempting an 
answer, it is helpful define the critical point with somewhat greater precision. In the PCP 



the bulk critical point p c is defined such that in the infinite-size limit, the stationary order 
parameter (density of NN pairs) is zero for p > p c (subcritical regime) and scales as p ~ A' 3 
for A = p c — p ~ 0. Now consider a spreading process, in which activity is initially localized 
(for example, at a single NN pair in the PCP). The subcritical regime for spreading is defined 
as the set of p values for which the survival probability P(t) — > as t — > oo; the critical 
value for spreading is then 

p c (0) = min{p : lim Pit) = 0} . (2) 

t— >oo 

In simulations, p c ((fi) is identified as the value for which the survival probability, mean 
population n(t), and mean-square spread R 2 (t) all follow asymptotic power laws, which 
define the exponents 5, rj and z. (In the PCP, n(t) and R 2 refer to NN pairs; in the 
simpler CP, of course, they represent the number and spread of particles. Note that in the 
supercritical regime, n ~ t d and R 2 ~ t 2 , that is, the finite fraction of trials that survive 
spread with a finite velocity, and have a finite bulk density For A ~ 0, we expect a crossover 
from the exponents r\ and z at short times to the supercritical growth laws at longer times.) 
We are not aware of any proof that P, n, and z follow power laws at p c {<P)- But it is known 
that in the CP, away from the critical point, the approach to the stationary state, be it 
active or absorbing, is exponential [37,41]; in some simpler cases, such as compact directed 
percolation, a power-law can be demonstrated explicitly [42,43]. (On the other hand, the 
CP with quenched disorder appears to have a logarithmic time-dependence at the critical 
point, and power- law away from the critical point [44,45]). In any event, we shall assume 
that in the PCP, p c {4>) can be identified via the power-law criterion. 

Consider the mean population n s (t), in trials that survive until (at least) time t. By 
definition, n s (t) > 1, for any value of p, but for p > p c (0), n s (t) remains bounded as t — > oo, 
while for p < p c (4>) it grows without bound. Unbounded growth also obtains at p c {4>) if we 
assume that P(t) decays slower than exponentially at the critical point (absenting a proof, 
we merely observe that exponential decay would be incompatible with scale-invariance!). 
That is because the rate of extinction is nronortional to the nrobabilitv of having exactlv 



one pair: 

~=pP(t)Pr[ n< (t) = l]. (3) 

If n s (t) remains bounded, then Pr [n s (t) = 1] > 0, and Eq. 3 implies exponential decay 
of the survival probability. Assuming the power laws noted above for P and n, we have 
ns{t) ~ t v+s ; unbounded growth of n s (t) at p c (4>) implies rj + 5 > 0. 

The above observations are helpful in analyzing the implications of a critical point shift. 
Suppose that p c (<f>) < Pc for spreading in an environment with < nat . (That is, to 
compensate for the "hostile" environment, the spreading critical point shifts to a value 
that lies in the bulk supercritical regime.) This means that for p c (<j>) < P < p c , the bulk is 
supercritical, but the spreading process always dies. At the critical point p c (4>), the spreading 
trials that do survive for long times have a nonzero bulk density p& ~ [p c — p c ((j>)]^, since 
the interior of the active region must eventually attain the stationary state. At long times, 
this interior region, having a finite density, makes the dominant contribution to n s and R 2 , 
implying that n s ~ (R 2 ) d / 2 ~ t zd ^ 2 , which in turn implies the hyperscaling relation 

V + 5 = — , (4) 

rather than Eq. (1). [Eq. (4) describes spreading at first-order transitions between different 
absorbing states, for example in compact DP or the voter model [46]. In this case Eq. (4) 
follows directly from Eq. (1), since (3 = 0.] 

Suppose, conversely, that p c ((f>) > Pc for > (f) na t, i-e., the "friendly" environment 
permits a spreading trial to survive indefinitely even if the bulk is subcritical. If so, then for 
Pc < P < Pc{4>), a surviving trial must have its activity concentrated in an annular region (or 
at two fronts, in one dimension), since the order parameter vanishes in the bulk. (One may 
think of the activity as a "chemical wave" that converts one kind of absorbing configuration, 
with > (pnat, into another, with ~ <f> na t-) In one dimension, moreover, the width of this 
active zone must grow fast enough that n s (t) grows without bound as t — > oo. 

Summarizing, if the critical point for spreading depends on the surrounding particle 
densitv. then we should exnect. for ch < ch„„t. that the activitv nrofile n„(x.t) in surviving' 



trials is plateaulike, and that the spreading exponents satisfy Eq. (4); for <fi > (p nat , we 
would expect the activity profile to be bimodal, with an inactive central region. There are 
no arguments, to our knowledge, prohibiting such behaviors. (Indeed, such modifications of 
usual spreading were reported for two-dimensional models in Refs. [28] and [29].) But if the 
activity profiles do not take the anticipated forms, it is likely that either no shift in p c (4>) 
exists, or that the asymptotic long-time behavior has not been probed. 

IV. SIMULATION METHOD 

Initial configuration. For each trial, we construct an initial configuration with a single 
NN pair or "seed" at the center of the system (sites L/2 and L/2 + 1 on a line of L sites). We 
consider three values of 0: zero, 1/2, and the natural density 0.242. For = 1/2 the sites 
immediately to either side of the seed are vacant, beginning an alternating sequence that 
frames seed. For the studies at <p nat we generate random initial configurations as follows. 
Consider the sites i > L/2 + 1. Since L/2 + 1 is occupied, it must be followed by a gap of 
g > 1 vacant sites before the next occupied site. We generate a sequence of gaps gi — l + m, 
where the are a set of independent Poisson random variables with mean (m,) = 4>nlt ~ 2, 
in order to reproduce the natural particle density. (The same procedure is applied to sites 
L, L — 1, ...1. In each trial we keep a record of the extent of the active region. To prepare 
for the following trial we only need to reset sites within that region; since it is usually much 
smaller than the full system, this results in a considerable speedup.) 

Dynamics. We maintain a list of the N p current nearest-neighbor pairs. At each step 
we choose a pair at random from the list, and a process (annihilation with probability p, 
creation with probability 1 — p). In case of annihilation, the pair is simply removed. For 
creation, we choose a site % at random from the two sites neighboring of the pair, and place 
a new particle at i if it is currently vacant. (If i is occupied the configuration remains the 
same.) The time increment associated with a step is At = 1/N P , corresponding to one 
transition per pair per unit time, in agreement with the transition rates that define the 



process. Following each change we update the list of pairs. The trial continues until either 
no pairs remain (extinction) or a preset maximum time t m is reached. (In the largest studies 
reported here, t m = 22026 ~ e 10 .) We use a lattice of L = 5 x 10 4 sites, sufficiently large 
that activity never reaches the boundary of the system for t <t m . 

Reweighting scheme. We have realized a significant speedup, permitting larger samples 
and longer durations than in previous studies, by means of a reweighting technique. The 
basis and details of this method are given in Ref. [36]. Briefly, we run the simulations at a 
central value, p , of the annihilation probability p, and reweigh the sample to study nearby 
values p' = p + nAp. (We always use p very near p C) typically p = 0.07709 or 0.0771; Ap 
is either 10~ 5 or 5 x 10~ 6 , depending on t rn .) Suppose that in a particular trial, there have 
been n a annihilation events and n c creation events (successful or not) up to time t. The 
weight for this sequence of events in a simulation with annihilation probability p', is 



which represents the ratio of the probability for this sequence in the process with parameter 
p' to the corresponding probability for p . In practice, we record the number of pairs, 
mean-square displacement, and continued survival of a trial at unit time intervals. The 
simulation estimate for the central value, po, of any property A at time t is, as usual, A t . po = 
N^. 1 J2k=i Ak,u where A^t is the value at time t in the kih of a total of N tr independent 
trials. The corresponding estimate for parameter value p' is A t . p > = N^ 1 J2k=i ^(p'iPo', t)A^t- 
The reweighting factor uj varies from trial to trial, depending on the number of creation and 
annihilation events realized up to time t. An analysis of the range of p values for which 
reweighting is useful is given in Ref. [36]. 

Application of this method to the CP yielded results of unprecedented precision for 
critical point and exponents [36]. Beyond the obvious economy of studying various parameter 
values at once, the fact that the results for all p values come from the same data set eliminates 
run-to-run fluctuations that tend to obscure the value of p c . (Our method, of course, does not 
eliminate fluctuations ner se. But here the fluctuations affect the results for all narameter 




(5) 



values uniformly. Without reweighting, one has to generate a different sample, with its 
particular fluctuations, for each parameter value of interest.) 

V. RESULTS 
A. Critical Point and Spreading Exponents 

We determine the survival probability, mean number of pairs, and mean-square distance 
of pairs from the original seed, expecting that these follow asymptotic power laws at the 
critical point, and show deviations from power laws for off-critical values of p. To identify 
Pc{4>) we follow the widely-used practice of studying the local slopes, S(t), rj(t) and z(t), 
defined as the derivatives of InP, Inn, and InP 2 , respectively, with respect to Int. We 
evaluate S(t) by performing a least-squares linear fit to the InP data for a set of 25 equally- 
spaced values (an increment of 0.1) of lnt; i](t) and z{t) are obtained similarly. 

In studies of spreading at an absorbing-state phase transition, the local slopes are usually 
plotted versus t~ l . The curves for 5(t) and r](t) typically fall into three groups: those that 
curve sharply upward are taken as marking supercritical parameter values, those curving 
sharply downward are associated with the subcritical regime, and those that seem to have 
a finite limiting value (as t^ 1 — > 0) are consistent with criticality. (The extrapolated value 
furnishes an estimate for the associated exponent, 5, i], or z. There is often some uncertainty, 
i.e., two or more parameter values may "look" critical, especially when fluctuations are 
strong. The curves for z(t) generally provide much less information about the location of 
the critical point.) The rationale for this practice is that in the one-dimensional CP the 
leading correction to scaling appears to be oc that is, for f > 1, 

5(t) - 5[1 + ar Al + bt~ A2 + •••], (6) 

with Ai = 1. In other cases, however, the dominant correction to scaling may decay more 
slowly than 1/t [47]. 



In Fig. 1 we plot the local slopes 5(t) from our spreading studies with = 0. (The data 
are from a set of 2 x 10 6 trials, extending to t m = 22026. The simulations were performed 
at p = 0.07709, and reweighted to study a series of p values at intervals of Ap = 5 x 10~ 6 .) 
When plotted versus as in Fig. 1, all of the curves swerve sharply downward; there is 
no obvious candidate for p c . At this point we recall an important observation of Odor et al.: 
the particle density <fi(t) in the active region approaches its natural value as a power law, 
\4>{t) — 4> na t\ ~ t~ 5op , with the standard DP exponent 5dp — 0.16 [35]. Given the coupling 
between the order parameter and the particle density, it is natural to suppose that the slow 
relaxation of the latter will yield a dominant correction to scaling ~ t~ &DP rather that ~ t~ x . 
We test this hypothesis in Fig. 2, where the 8{t) data are replotted versus t~ 5op . Evidently, 
the strong curvature of the 5(t) plot in Fig. 1 is accounted for by a change in the correction 
to scaling exponent; the graphs of Fig. 2 are roughly linear. Similar changes attend the 
switch from t^ 1 to t~ Sop in the other local-slope plots, for all of the values studied here. 

In Fig. 2, we see that the curve for p = 0.07707 is nearly linear, whilst those for 0.077065 
and 0.077075 show noticeable curvature. (This visual impression is confirmed by least- 
squares quadratic fits to the data.) Linear extrapolations of the 5(t) plots yield 5 = 0.252, 
0.256, and 0.261 for p = 0.077065, 0.07707, and 0.077075, respectively. Based on these 
results, we estimate 5^ =0 = 0.256(4). Fig. 3 shows a similar plot for r](t), the local slope of 
the population size. The curves are somewhat less regular, and the most linear plot appears 
to be that for p = 0.077075. Linear extrapolations for p = 0.07707, 0.077075, and 0.07708 
yield rj = 0.218, 0.211, and 0.204, leading to the estimate r^ =0 = 0.211(7). The data for z{t) 
are plotted in Fig. 4. Here it is difficult to discern p c ; linear extrapolations over the range 
p = 0.077065 - 0.07708 yield 2^ =0 = 1.245(7). Our results are consistent with those reported 
in Ref. [25] for spreading in an initially empty lattice: <5^ =0 = 0.250(5), r/^Q = 0.215(5), 
and z^=o = 1.238(4). (The somewhat larger error bars in the present results reflect the 
uncertainty inherent in extrapolating to t~ - 16 = 0.) Combining the critical point estimates 
from the analyses of 5(t) and r](t), we obtain p c (0) = 0.077073(6). 

Next we consider snreadinff in a, half-filled environment. Here we nerformed a, total of 



1.2 x 10 6 trials, extending to t m = 22026; in this case we used p = 0.0771 as the central 
simulation value. The local slopes 5(t) and rj(t) are shown in Figs. 5 and 6, respectively. 
Note that —5 and rj approach their limiting values from below, while for = the approach 
is from above. Analyzing the local slopes as above, we find p c (l/2) = 0.077087(6). The 
exponent estimates are listed in Table I. Finally, we performed simulations of spreading into 
an environment with particle density 4> nat = 0.242. The central value of p was again 0.07710, 
and the number of trials 2 x 10 6 , but the studies extended only to time 8103. Analyses of6(t) 
and rf(t) (the latter is shown in Fig. 7) yield p c {4>nat) = 0.077093(7). The critical exponents 
for = (finat, listed in Table I, are in good agreement with the known DP values. 



B. Activity and particle density profiles 

The discussion of Sec. Ill motivates our study of the spatial distribution of activity (the 
pair density) and the particle density in surviving trials. According to the scaling hypothesis 
[34], at the critical point, the mean order-parameter density p s (x, t) in surviving trials, given 
an initial seed at x = and t — 0, should scale as 

p a (x,t)=f' 2 - s -'> R(x 2 /t z ), (7) 

where the scaling function R vanishes for large values of its argument. The hyperscaling 
relation Eq. (1) implies that the prefactor is t~ &DP . We therefore expect plots of p*(x*,t) = 
t Sop p(x/t z l 2 , t) to collapse onto a single curve. Figures 8, 9, and 10 show such plots for = 0, 
1/2, and 4> nat = 0.242, respectively, confirming the anticipated scaling. (The study for = 
was performed using p = 0.07708, and those for = 0.242 and 1/2 using p = 0.07709; there 
is no rewighting.) The scaling function has the same qualitative form in all three cases; 
there is no hint of the plateaulike or bimodal profiles envisaged in Sec. III. The fact that 
the scaled activity profile is stable for times > 1000 suggests that the local slopes 5(t), rj(t) 
and z(t) have also attained their asymptotic scaling regime, so that the critical exponents 
can be extrapolated with confidence. 



We note also that the rate of spreading is greater, the larger is 0. This can be quantified 
by measuring the mean-square spread; we find x* 2 ~ 0.75i, I. It, and 2.0t for = 0, 0.242, 
and 0.5, respectively. Most significantly, the definition x* = xjt z l 2 involves the 0- dependent 
z values listed in Table I. Were we to use, for example, the standard DP value z = 1.265 in 
x* for = we would find that x* 2 /t approaches zero as t — > oo. 

A similar analysis can be applied to the mean particle density (f> s (x, t) in surviving trials. 
Here 4> s (x,t) is fixed at its initial value for x t z ^ 2 , that is, outside the active region, 
and we expect <p s {x,t) — 4> na t to decay ~ t~ Sr,p within the active region. These behaviors 
are reflected in the scaling form 

+ r SDP F{x 2 /t z ) . (8) 

Here the first term on the r.h.s. represents the uniform particle density outside the active 
region; the factor in square brackets switches from 1 outside, to zero well inside this region. 
(For = (fi n at, of course, this term vanishes.) The second term represents relaxation of the 
particle density to its natural value within the active region. To test this scaling hypothesis 
we plot 

<f>*(x*,t) EE t SDP |(0 s (x,t) - <j) nat ) - (0 - nat ) 

versus x*. Figures 11, 12, and 13 are scaling plots of (p*(x*,t) for = 0, 1/2, and 0.242, 
respectively. The data collapse confirms the observation of Odor et al. that the particle 
density exhibits a power-law relaxation to its natural value, |0(i) —(p na t ~ t~ &DP \, within the 
active region [35]. 

VI. DISCUSSION 

We return to the three clues mentioned in the Introduction. First our simulations, which 
are the most extensive (in terms of sample size and trial duration), of the PCP to date, 
shows no evidence for a shift of the critical point with particle density. It is true (see Table 



4>s (X, t) - nat = (0 - (j) na t) 



Ps(x,t) 
Ps(P,t) 



Ps(x,t) 
Ps(0,t) 



(9) 



I) that our p c estimate for = 1/2 is slightly higher than that for = 0, but the data for the 
three values studied show no systematic trend. In particular, our data do not confirm the 
results reported in Ref. [35], where p c is given as 0.07704 for = and 0.07714 for = 0.432 
(and would presumably be at least as large, for = 1/2), so that p c (l/2) — p c (0) > 10~ 4 . 
Here, by contrast, we find p c (l/2) -p c (0) = 1.4 x 10~ 5 , with an uncertainty of ±1.2 x 10~ 5 . 
To summarize, the present results are incompatible with the relatively large shift reported 
in Ref. [35], but are quite compatible with there being no shift at all. 

In Figures 14, 15 and 16 we plot our results for the spreading exponents 5, i], and z as 
functions of 0, together with results from Refs. [25] and [35]. The overall consistency among 
these studies leaves little doubt that there is a systematic, roughly linear dependence of the 
exponents on the particle density. 

Our results for the spreading exponents are sufficiently precise to rule out the proposal 
that z and 5 + 77 are independent of 0. While their variation is much smaller than that in 
5 or r] separately, there are significant changes, amounting to a relative variation of about 
5% between the extremes = and = 1/2. The fact that z varies with was confirmed 
in our scaling analysis of the activity profile. Finally, as shown in Table I, the generalized 
hyperscaling relation is satisfied to within numerical uncertainty. Given the variation in z, 
a compensating variation in 5 + i] is needed for hyperscaling to hold. That 5 + r\ and z 
depend on forces us to abandon what seemed an attractive simplifying hypothesis: the 
environment affects not only the likelihood of survival, but (to a lesser extent) the scaling 
properties of trials that do survive. 

Can we understand nonuniversal spreading in models with INAC? Thus far, attempts 
to identify a marginal parameter in a renormalization group analysis of the PCP and allied 
models have not borne fruit. Nonuniversality does appear to be associated with a long 
memory, due to the slow relaxation of an auxiliary field coupled to the order parameter. Such 
a memory term appears in a field theory that successfully predicts the static critical behavior 
of the PCP [27], and that reproduces (numerically) nonuniversal spreading [32]. Moreover, 
varvinff exponents have been observed in one-dimensional models with a, nninne absorbing 



configuration and a slowly-relaxing memory of the initial condition [29,30]. Explaining the 
nonuniversality of the spreading exponents, and predicting their values as a function of the 
particle density, remain as theoretical challenges. 
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TABLES 

TABLE I. Critical point p c (4>) and spreading exponents for the PCP. DP exponents from Ref. 
[48]. Numbers in parentheses denote uncertainties in the last figure. 



4> 


Pc(4>) 


5 


V 


z 


5 + 7] 


2 V + 2(5 + 5 DP ) - z 





0.077073(6) 


0.256(4) 


0.211(7) 


1.245(7) 


0.467(11) 


0.01(3) 


0.242 


0.077093(7) 


0.161(2) 


0.314(4) 


1.264(4) 


0.475(6) 


0.005(16) 


0.5 


0.077087(6) 


0.056(3) 


0.443(5) 


1.303(3) 


0.499(8) 


0.014(19) 


DP 




0.1595 


0.3137 


1.2652 


0.4732(1) 





Figure Captions 

Fig. 1. Local slope —S(t) versus t^ 1 for the PCP with = 0. The central curve (with data 
points) is the simulation result at p — 0.07709; curves above and below were obtained via 
reweighting, for intervals Ap = 5 x 10~ 6 . The inset is a detail of the late-time behavior. 

Fig. 2. The data of Fig. 1 plotted versus t~ 016 . 

Fig. 3. Local slope rj(t) for = 0, versus t~ 016 , curves as in Fig. I. 

Fig. 4. Local slope z(t) for = 0, versus t~ 016 , curves as in Fig. I. 

Fig. 5. Local slope —5(t) versus t -016 , for = 1/2. The central curve (with data points) is 
the simulation result at p — 0.07710; curves above and below were obtained via reweighting, 
for intervals Ap = 5 x 10~ 6 . 

Fig. 6. Local slope rj(t) versus t -016 , for = 1/2; curves as in Fig. 5. 

Fig. 7. Local slope rj(t) versus t~ 016 , for = 0.242. The central curve (with data points) is 
the simulation result at p — 0.07710; curves above and below were obtained via reweighting, 
for intervals Ap = 5 x 10~ 6 . 

Fig. 8. Scaled activity density p* = t s ° p p versus x* = xjt z l 2 in the critical PCP with = 0. 
Open squares: t = 1000; diamonds: t = 5000; grey squares: t = 10000; line: t = 20000. 

Fig. 9. Scaled activity density as in Fig. 8, but for = 1/2. 

Fig. 10. Scaled activity density as in Fig. 8, but for = 0.242. 

Fig. 11. Scaled particle density 0* versus x* in the critical PCP with = 0. Symbols as in 



Fig 



8 



Fig 



12 



Scaled particle density as in Fig. 11, for = 1/2. 



Fig. 13. Scaled particle density as in Fig. 11, for = 0.242. 

Fig. 14. Critical exponent 5 versus particle density 0. Filled squares: present study; open 
squares: Ref. [35]; circles: Ref. [25]. 

Fig. 15. Critical exponent rj versus particle density 0. Symbols as in Fig. 14. 
Fig. 16. Critical exponent z versus particle density 0. Symbols as in Fig. 14. 
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